Quintom cosmologies admitting either tracking or phantom attractors 
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In this paper we investigate the evolution of a class of cosmologies fuelled by quintom dark energy 
and dark matter. Quintom dark energy is a hybrid of quintessence and phantom which involves 
the participation of two reals scalar fields playing the roles of those two types of dark energy. In 
that framework we examine from a dynamical systems perspective the possibility that those fields 
are coupled among them by considering an exponential potential with an interesting functional 
dependence similar but not identical to others studied before. The model we consider represents a 
counterexample to the typical behavior of quintom models with exponential potentials because it 
admits either tracking attractors (w = 0), or phantom attractors (w < — 1). 

PACS numbers: 98.80.-k, 95.36.+X, 98.80.Jk 



I. INTRODUCTION 

Understanding the nature of dark energy has become a 
joint non trivial pursuit of an everyday growing number 
of cosmologists. It is expected that the large investment 
which is being made in technical and human resources 
will pay off in the form of a deeper understanding of the 
link between the physics governing the universe at the 
classical and quantum levels. 

Abundant observational evidence supports the view 
that the universe is currently undergoing accelerated ex- 
pansion 0, 0, 0, 0, HI El due to a time-evolving compo- 
nent of unknown nature called dark energy. Peculiarly, 
in most of the analyses the best agreement with obser- 
vational data is provided by models in which, as time 
goes forward, the dark energy equation of state param- 
eter w = p/p crosses the w = — 1 divide (to become 
even more negative). However, among recent studies us- 
ing newest SN data, some suggest the convenience of the 
phantom divide crossing may just be an illusion due to 
systematic errors in observations Q, whereas some oth- 
ers conclude that all classes of dark energy models are 
comfortably allowed 8]. The number of studies with a 
negative regard on the existence of aw < —1 epoch are 
still scarce, and to add to the problem results have strong 
dependence on the chosen parametrization of w in terms 
of redshift |!j , so there is not a full consensus about the 
observational preference of the phantom divide crossing. 
For this reasons, we feel much has to be done yet before 
one can definitely say the scale tilted to one side, so con- 
tinuing to get deeper insight into models which do the 
crossing seems still justified. 

Some interesting general aspects of the pro blem of the 
phantom divide crossing were discussed in [l(J] , where the 
viability requirements on the equation of state and sound 
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speed were analyzed. Even though, some of realizations 
of the crossing can have an extradimensional origin, ei- 
ther in the brane or the string gas context 0] , scalar 
field models in standard four dimensional physics are the 
most popular options of the inventory. The impossibility 
of the occurrence of the transition in traditional single 
field models 0| has motivated much activity in the con- 
struction of two field models that do the job. Examples 
of explicit constructions can be found in [hi [Til l2l| , but 
perhaps the class of models which have received most at- 
tention are quintom cosmologies 0, 0, 0, 0, |2(], |2l| . 

Other worth mentioning alternatives are models in the 
framework of scalar-tensor theories , a single field pro- 
posal involving high order derivative operators in the la- 
grangian |23| |. a model with a single dynamical scalar 
field coupled to an a priori non-dynamical background 
covariant vector field [24J, and an interacting Chaplygin 
gas [25I . We stick here to the main stream and consider 
quintom cosmologies as well. 

The cosmological evolution of quintom models has 
been examined using standard dynamical systems tech- 
niques in ^1 and pjjj. In these references an isotropic 
and homogeneous universe was supposed to be filled with 
dark matter and dark energy, and the potential of the 
latter was chosen to be of exponential form. The differ- 
ence between and |20l | is that in the second case the 
potential accounts for an interaction between the con- 
ventional and the phantom scalar fields. Here we revisit 
quintom models from the dynamical systems perspective, 
but make a choice of potential which is closely related to 
that in |20| . but which at the same time does not in- 
clude it as a particular case. The subtle difference of our 
choice has unexpected consequences. Firstly, there exists 
the possibility that quintom dark energy (with fairly gen- 
eral initial conditions) tracks dark matter, thus meaning 
that avoidance of tracking behavior in quintom models 
is not generic. Secondly, in some cases for which the pa- 
rameters in the potential exclude the tracking attractor, 
then the model may have a purely phantom attractor, 
i.e. this will be models in which w relaxes to be less 



2 



than — 1. This is the converse of the behavior observed 
in previously studied models with exponential potentials, 
in which the late-time asymptotic behavior corresponded 
simply to w = — 1, despite the presence of a transitory 
epoch with w < — 1. In our opinion, this finding shows 
the casuistics of the evolutionary behavior of quintom 
dark energy models is in fact reacher than was known 
up to the date, and makes it a subject of study worth of 
further investigations. 

II. THE MODEL 

We investigate the evolution of a spatially flat 
Friedmann- Robertson- Walker (FRW) spacetime filled 
with dark matter with energy density p m and quintom 
dark energy with energy density p de and pressure p de . 

Our model is assumed to obey the standard Friedmann 
equation and energy conservation equations, that is, 

3H 2 = p rn + p de , (1) 
p de + 3H(p de + Pde ) = 0, (2) 
p m + 3Hp m = . (3) 

Here and throughout overdots denote differentiation with 
respect to cosmic time t, H = a/ a is the Hubble factor, 
and a is the synchronous scale factor. 

Combining Eqs. 11121 |3J one obtains the evolution 
equation of H, i.e., 

-2H = p m + p de + p de . (4) 

Dynamical systems are a well-known tool for depicting 
the asymptotic (and sometimes the intermediate) behav- 
ior of cosmological models by making use of the concepts 
of past and future attractors. We, thus, take advantage 
of that fact and apply it to the model under considera- 
tion. One of the first tasks involved in this protocol is 
introducing a set of convenient (expansion normalized) 
variables which allow casting the conservation equations 
and the evolution equation of H as a dynamical system, 
subject to a constraint arising from the Friedmann equa- 
tion 0J. 

In our quintom cosmologies a minimum of four vari- 
ables are needed to construct an autonomous dynamical 
system, but because of the existence of the constraint 
one can always "forget about" the evolution of one of the 
variables. Let us consider just for a while we are using n 
variables which we will denote as x' . The equations will 
be of the form 

x*' = f{x\x\...,x n ). (5) 

The prime are used to denote differentiation with respect 
to an alternative new time variable r which has to be 
chosen upon convenience dictated by the normalization 
of the variables (here a very standard choice will suffice). 

The next step in the study of the evolution of our dy- 
namical system is to find its fixed (or critical) points 



(x 1 *, x 2 *, . . . , x™*), which are given by the conditions 

f(x i *,x 2 *,...,x n *) = a. (6) 

The stability of the fixed points (x 1 * , x 2 * , . . . , x n *) is then 
analyzed by studying the linearized dynamical system 
obtained by expanding the n evolution equations about 
those fixed points (as explained in many seminal refer- 
ences, e.g 26j). Afterward, one tries solutions in the 
form (x 1 , x 2 , . . . , x n ) = (ci, C2, . . . , c„) e At in the linear 
approximation, and finds that their characteristic expo- 
nent A and the constant vector (ci, C2, . . . , c n ) must be 
respectively an eigenvalue and an eigenvector of the ma- 
trix 
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(7) 

Let us know describe the basics of the character of the 
fixed points depending on the values of the characteris- 
tic exponents. If the real part of the three characteristic 
exponents is negative, the fixed point is asymptotically 
stable, i.e., an attractor. On the other hand, it is enough 
to have (at least) one characteristic exponent with pos- 
itive real part to make the fixed point asymptotically 
unstable, i.e., a repeller. This repeller is a saddle point 
if at least one of the other characteristic exponents has 
a negative real part, in which case there is, apart from 
the unstable manifold, a stable manifold containing the 
exceptional orbits that converge to the fixed point. The 
case in which the largest real part is precisely zero must 
be analyzed using other methods, because the linear anal- 
ysis is inconclusive. The geometric form of the orbits near 
the fixed point is determined by the imaginary part of the 
characteristic exponents. If the three are real the fixed 
point is a node. A couple of complex conjugate exponents 
leads, except in degenerate cases, to a spiral saddle (the 
orbits are helices near the fixed point). In addition, when 
one of the exponents is null the point is not hyperbolic 
and therefore structural stability cannot be guaranteed 
(the geometric form of the trajectories may change under 
small perturbations). 

In the next sections we are going to present our quin- 
tom model and will study its evolution equations. We 
will find and characterize the fixed points and we will 
analyze the cosmological consequences of the results ob- 
tained. Special attention will be paid to the possibility 
these models offer to explain the passage from a conven- 
tional to a phantom equation of state. 
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TABLE I: Location, existence and decelation factor of the critical points for m > 0, n > and y > 0. 



Name x& y Existence 



O All fh and n 



C± ±<Jl + x* 2 x* v All fh and n 2 



P fh -n ^l-fh 2 + h 2 m 2 -n 2 <l -l + 3(m 2 -n 2 ) 

m — n > 1/2 



2(m 2 — n 2 ) 2(m 2 — n 2 ) 2\/m 2 — n 2 2 



III. PHASE-SPACE 

We begin by presenting Eqs. ill'IMI II) in the form of a 
dynamical system. This is done by making the following 
convenient choice of variables: 

- t - t - - V^" (o\ 

which renders the Friedmann equation as 

x\-x% + y 2 + z 2 = l. (9) 

Eq. m leads to the constraint @ which allows con- 
sidering only the evolution of those three variables solely 
(the evolution of the fourth one been determined by the 
evolution of the former). In this case we will choose z as 
the variable one can mostly do without. 

The choice of potential is crucial to our discussion, we 
will set 

7 = y e-^^+ ft ^, (10) 

with fh, h and Vq positive constants. This kind of po- 
tential was considered to account for the interaction of 
conventional scalar fields for instance in p^|. 

In this paper we will only be concerned with expand- 
ing universes, that is, in what follows we will restrict 
ourselves to H > or equivalently to y > 0. Combining 
expressions (llll'l.'il 1181) . the following evolution equations 
are obtained: 

a£ = i (3T™, 2 + (q - 2)x 4 ) (11) 

x; = -~(3rm/ 2 -(<Z-2K) (12) 

y' = ^(1 + q - 3(mx,p + hx v ))y (13) 

Here primes denote differentiation with respect to a 
new time variable r = log a 3 and q = —da/ a 2 stands for 
the deceleration factor. Explicitly 

q = \ (3 [x 2 - x% - y 2 ) + l) . (14) 




FIG. 1: Phase-space trajectories of the 3D system for m = 0.5 
and n — 0.6. 




FIG. 2: Phase-space trajectories of the 3D system for m = 
0.75 and n = 0.05. 



The evolution equations of variables x^ , x v and y form 
a 3D dynamical system defined on the phase space 

* = {(x4,x v ,y) :0<xl-xl+y 2 < 1}. (15) 

Thus, the phase space trajectories live in an hyperboloid. 

The linear analysis described above tells us the 3D dy- 
namical system under consideration has at most three 
isolated points (depending on the values of fh and h). 
We will denote them as O, T and P. In addition, there 
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are two curves (hyperbolas) of non-isolated fixed points, 
and we will denote them as C±. 

Some information regarding the dynamical character 
of those fixed points is presented in Tables i| I II II) . Be- 
low we will complete that information with an identifica- 
tion of the cosmological models represented by the fixed 
points of our system, and in the case of the isolated fixed 
points we will also comment on their geometry. Finally, 
we will support our findings with numerical simulations. 
Note that in the case of 3D dynamical systems it is in 
general difficult to retrieve enough information on the 
system without resorting to numerical experiments |27| . 
The need to employ numerical tools is particularly man- 
ifest in the problem we are concerned with, because the 
asymptotic and intermediate behavior stand on the same 
basis of importance. Bear in mind that we do not only 
want to show our system admits solutions with equation 
of state parameter w > — 1 and other with w < — 1, but 
we also want to show that a transition between such so- 
lutions can occur without too much fine-tuning in initial 
conditions. In our case 



w 



V 



K + v 



(16) 



and we will evaluate the latter at the fixed points. 

In broad terms, if one of the asymptotic limits of 
our model is to give good description of the universe at 
present it must be an accelerating solution. Thus, non- 
accelerated solutions must be asymptotically unstable, or 
put another way, they must be disfavored by initial con- 
ditions. In this sense only the cases in which P is stable 
would be satisfactory. 

As can be deduced from Tables dill II) . many of the 
features of our dynamical system depend on the value of 



the quantity m 



One can easily realize why such 



quantity may be of relevance just by noticing 



6(m 



V 



2 /T/ \ 2 

V 



(17) 



i.e. the quantity marking the existence of one attractor 
or the other compares the slope of the potential in two 
different directions, or equivalently it compares how fast 
the fields release (gain) potential energy when they roll 
down (climb up) the potential. Unlike in single scalar 
field models with exponential potentials, one can have 
accelerated expansion even if the potential is not too shal- 
low, what matters here is not how flat the potential is, 
but rather whether it is much flatter in one direction that 
in the other. 

After these preliminary remarks, let us now be more 
specific about the dynamical character of the isolated and 
non-isolated fixed points of our system. 

The point O represents a matter dominated deceler- 
ating solution. It is a saddle, and its unstable character 
(see all four figures) agrees with what one might have an- 
ticipated. Such solutions are only expected to be relevant 
at early times. The variables we are using do not allow 




FIG. 3: Phase-space trajectories of the 3D system for m = 0.9 
and n — 0.4. 




FIG. 4: Phase-space trajectories of the 3D system for m = 2 
and n = 0.5. 



to determine the value of equation of state parameter 
corresponding in this case (because we do not know the 
values of the different ratios between x^, x v and y). This 
task would possibly require defining an alternative set of 
coordinates more suitable for an early time description, 
but this is of little interest in the context of this paper. 

The curves C± represents solutions in which the con- 
tribution of matter and the potential energy to the total 
energy density is negligible. These solutions are therefore 
of stiff-fluid type (w = 1), which in turn correspond to 
a decelerating universe. The unstable character of this 
solutions means they are disfavored from the initial con- 
ditions point of view, and therefore it is unlikely they can 
represent the final stage in the evolution of our universe. 

Further interesting remarks concerning the curves C± 
are in order to help anticipate the intermediate behavior 
of our dynamical system. In the invariant set y = which 
contains the point O as well as the curve C± we have 
identified the monotonic function /C, which is defined as 



K 



(18) 
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TABLE II: Eigenvalues, dynamical character and equation of state parameter of the fixed points for m > and n > 0. We use 
A = ^/-7 + 4/(m 2 - W) . 



Name Eigenvalues Dynamical character 

O ~ 5' 5 unstable undefined 



C± ( 1, 0, 1 — ni* =p roJl + a;* 2 I unstable 1 



P (-1 + 2(m 2 — fi 2 ), — 1 + m 2 — n 2 , — 1 + m 2 — n 2 ) stable if m < J i + rt 2 , unstable otherwise -1 + 2(m 2 

-i,- i(l + A), - A) J stable if either A 2 < or if 1 > A, 

unstable otherwise 



It satisfies 

K! = -2K, (19) 

so as expansion proceeds K, decreases. This monotonic 
function takes its maximum and minimum values at x 2 ^ — 
XT, = 1 and Xcf, — 0, so typically trajectories will transit 
between C± and O, and one will be able to say the curves 
C± are local sources. This appreciation is supported by 
the four figures. 

The point T represents a solution in which quintom 
dark energy tracks matter (the equation of state of the 
quintom fluid is dust-like), and when such point exists 
it is an attractor. This late-time asymptotic state rep- 
resents a decelerating cosmological model in which the 
fractional energy densities of matter and dark energy are 
proportional. Given its eigenvalues structure this solu- 
tion is either a saddle or an spiral saddle. Since this 
solutions do not match the observed universe at present, 
a quintom model admitting such attractor is not satis- 
factory. Peculiarly, the existence of this solution is com- 
pletely inherent to the interaction between the two fields. 
There is no equivalent solution in the quintom models 
studied by. Its existence, in our setup at least, indicates 
quintom models admit other attractors than phantom or 
de Sitter ones, a fact which has not been noticed in the 
literature before (as far as we are concerned). 

The point P represents a solution in which quintom 
dark dominates over matter (the equation of state of the 
quintom fluid corresponds to a fluid that redshifts faster 
than dust), and when such point exists it is an attractor. 
This late-time asymptotic state does not necessarily rep- 
resents an accelerating cosmological model, that depends 
on the quantity fh 2 — n 2 . Given its eigenvalues structure 
this solution is either a saddle or a stable node. The ac- 
celerated solutions associated with this fixed point can 
provide a good representation of the observed universe 
at present. The potential we have chosen here allows 
this fixed point be characterize d by w < — 1, whereas in 
the quintom models studied by |19t |20| the attractor had 
W = — 1 necessarily. 

According to their eigenvalues, in principle the fixed 



points T and P could respectively have three and two 
different dynamical behaviors 1 . However, in some cases 
one point have a specific dynamical behavior enforces the 
nonexistence of the other point. Additionally, the exis- 
tence conditions of these points on one hand, and numer- 
ical analysis on the other, will help us identify the hetero- 
clinic sequences of these models. Under the assumptions 
m > and n > 0, four cases can be distinguished (see 
Fig. El): 

• Case i) For m < \Jfi 2 + 1/2, the point P is a stable 
node, whereas the point T does not exist. The het- 
eroclinic sequence in this case is C± — ► O — ► P 
(see Fig. 0). 

• Case ii) For ^Jn 2 + 1/2 < m < ^Jn 2 + 4/7, the 
point T is a stable node and the point P is a sad- 
dle. For these conditions the heteroclinic sequence 
is C± — ► O — ► T — ► P (see Fig. 0. 



• Case iii) For yjh 2 + 4/7 < rh < vl + n 2 , the point 
T is a spiral saddle and the point P is a saddle. 
For these conditions the heteroclinic sequence is the 
same as in the former case (see Fig. [3| . 

• Case iv) For m > yl + n 2 the point T is a spiral 
saddle whereas the point P does not exist. The 
heteroclinic sequence in this case is C± — ► O — ► 
T (see Fig. EJ). 

Before closing this section it is interesting to make a 
few remarks about the fractional energy densities Vt m = 
p m /3H 2 and fide = Pde/3H 2 - It is not difficult to see 
that 

f^n\' = w ^n (20) 



1 The point T can either be a saddle, a spiral saddle, or a stable 
node. The fixed point P can be a saddle or a stable node. 
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FIG. 5: Regions in the (n, m) plane representing the four 
cases of existence conditions and dynamical behavior of the 
fixed points T and P. 



This result can be viewed as a consistency proof of some 
of our previous results in the sense that in the course of 
the evolution dark energy becomes the dominant compo- 
nent when w < 0, because in this case it redshifts faster 
than dark matter. 




FIG. 6: Typical behavior of the fractional energy densities in 
the case where the phantom attractor exists. The dark and 
light gray lines depict fl m and Q,de respectively. 



Note as well that, unlike in the case of pure phantom 
dark energy, here the fractional energy densities of dark 
energy and dark matter remain both non-negligible for 
a long enough time span, as shown in Fig. |SJ So at 
least this problem which invalidates the simplest models 
with a late-time phantom behavior 29] is absent from 
our model. 

Needless to say that the trajectories we have presented 
show that for a wide range of initial conditions for which 
w ~ 1 (close to the hyperbolas), the region w < — 1 



is reached provided m 2 — n 2 < 0, so the crossing does 
indeed occur. 

IV. CONCLUSIONS AND FUTURE 
PROSPECTS 

The construction of possible realizations of a late-time 
crossing of the w = — 1 barrier in dark energy models 
has become an active area. Some of the studies take 
advantage of the dynamical systems approach to show 
such crossing is possible in a family of models dubbed 
quintom cosmologies. These are configurations with two 
scalar fields, and in these dynamical systems investiga- 
tions no other than exponential potentials have been cho- 
sen so far given some underlying physical motivation and 
the tractability of the choice. Here we stick to that 
simple election too, but consider a specification (still of 
exponential form) which is not a particular case of the 
choices made in other works. More specifically we have 
chosen a potential of the form V = V e -V6[m<f,+n V )/ 
with Vo, rh and n positive constants, whereas (follow- 
ing our notation for simplicity) related papers consid- 
ered either V = V 01 e - 2 ^* + y 02 e -2V6n*> []]| or y = 
Voi e - 2v ^™^ + V 02 e - 2VEnv + \V e -v^(™0+™ v ) 20], with 
Voii Vq2, and A positive constants as well. 

The potential we have chosen is responsible for an in- 
teraction between the fields and leads to some surprises. 
We have found that, contrary to previous beliefs, quin- 
tom cosmologies may admit tracking attractors. When 
the conditions for the existence of such attractors are 
not met, one has either phantom (w < — 1) or de Sitter 
(w = —1) attractors instead. This represents a departure 
from the situation described in 0] and j^, because in 
those cases the only attractors were de Sitter-like, al- 
though the trajectories ending there had w < — 1 tran- 
siently before entering the de Sitter phase. 

The model analyzed here can be criticized for instance 
on the grounds that our choice of potential is too spe- 
cific, but as happens for usual quintessence models, the 
results obtained here might have applications for related 
models with more general potentials. Another aspect of 
the problem is that the coincidence problem is not solved 
in this model; perhaps a way to address it is to consider 
interactions of the scalar fields with matter as well, but 
this is left for future projects. 
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